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two-dimensional square lattice 


Yong Zhu,! Ziqing Yang,'! Xin Zhang,' and Xiaosong Chen! 


‘State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 
Chinese Academy of Sciences, P.O. Box 2785, Beijing 100190, China 
(Dated: April 30, 2015) 


We have investigated both site and bond percolation on two dimensional lattice under the random 
rule and the product rule respectively. With the random rule, sites or bonds are added randomly 
into the lattice. From two candidates picked randomly, the site or bond with the smaller size product 
of two connected clusters is added when the product rule is taken. Not only the size of the largest 
cluster but also its size jump are studied to characterize the universality class of percolation. The 
finite-size scaling forms of giant cluster size and size jump are proposed and used to determine the 
critical exponents of percolation from Monte Carlo data. It is found that the critical exponents of 
both size and size jump in random site percolation are equal to that in random bond percolation. 
With the random rule, site and bond percolation belong to the same universality class. We obtain 
the critical exponents of the site percolation under the product rule, which are different from that 
of both random percolation and the bond percolation under the product rule. The universality class 
of site percolation differs different from that of bond percolation when the product rule is used. 
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I. INTRODUCTION 


Traditionally, percolation phase transitions are con- 
sidered to be continuous in various networks. How- 
ever, Achlioptas et al. concluded that the percola- 
tion transition in random network is discontinuous un- 
der the Achlioptas process (AP) in which the growth of 
large clusters is suppressed. During the evolution of a 
network under AP, edges are added into the network un- 
der product rule that the edge with minimum product of 
cluster sizes is connected from two randomly picked unoc- 
cupied edges. The percolation transitions under product 
rule in two-dimensional square lattice and scale- 
free networks are also investigated and argued to 
be discontinuous. Later, a series of investigations 
mi showed that different percolation transitions un- 
der product rule are actually continuous. In Ref. i, 
Costa et al. [7| studied a model with stronger suppression 
of large clusters and its percolation transition is contin- 
uous. Lee et al. concluded that the explosive perco- 
lation is continuous by studying cluster size distribution. 
Grassberger et al. [d] found that the explosive transition 
is still continuous but with unusual finite-size behavior. 
In the works of Liu et al. and Fan et al. H, a gen- 
eralized Achlioptas process (GAP) was introduced to in- 
vestigate the percolation transitions of two-dimensional 
lattice network and Erdös-Rényi (ER) network re- 
spectively. In the generalized Achlioptas process, the 
edge with minimum product of cluster sizes is connected 
with a probability p from two randomly picked unoccu- 
pied edges. When p = 1/2, GAP recovers to a random 
process. With p = 1, GAP becomes AP. The percola- 
tion transitions of two-dimensional lattice network and 
ER network under GAP keep to be continuous. Their 
critical exponents and and therefore universality classes 


depend on the probability parameter p. With a rigorous 
mathematical proof, Riordan and Warnke concluded 
that the percolation transitions under all Achlioptas pro- 
cesses are continuous. 


The site percolation on a two-dimensional square lat- 
tice under product rule was first investigated by Choi et 
al. mg. They claimed that the percolation is discontin- 
uous from the non-vanishing hysteresis between the di- 
rected and reverse process for the fraction of the sites be- 
longing to the largest cluster. However, Bastas et al. 
demonstrated that the hysteresis phenomena will disap- 
pear in the thermodynamic limit. Using the finite-size 
behaviors of the largest cluster and its standard devia- 
tion, they concluded that the explosive percolation tran- 
sition is continuous but belongs to different universality 
class. To clarify this controversy, further investigations 
are needed. 


The notation of a universality class is a basic tenet in 
the physics of critical phenomena. Within a universality 
class, the universal quantities such as critical exponents 
and scaling functions are independent of microscopic de- 
tails. It is accepted traditionally that a universality class 
is characterized by the spatial dimensionality of the sys- 
tem and by the number of the components of the order 
parameter. (See, e.g., the review article ie.) On a two- 
dimensional lattice, the random bond percolation has the 
same critical exponents as the random site percolation. 
The random bond and site percolation belong to the same 
universality class. For the bond percolation under GAP 
with p > 1/2, their critical exponents are different from 
that of random bond percolation and depend on the prob- 
ability parameter p . The universality classes of bond 
percolation under GAP are different from that of random 
bond percolation [o]. It is of interest to investigate the 
critical behaviors and universality class of site percola- 


tion in relation to the Achlioptas process. 

In this paper, we study the critical behaviors of site 
percolation under AP in two-dimensional square lattice. 
We investigate the sizes of giant clusters in the lattice. 
From the finite-size scaling behaviors of giant clusters, 
we can determine the percolation transition point and 
the corresponding critical exponents. Both the critical 
exponents of the sizes of giant clusters and their size gaps 
during evolution (3, [171 are calculated. 

Our paper is organized as following. In Section 2, we 
present the finite-size scaling behaviors of giant clusters 
near the phase transition point of percolation. The finite- 
size scaling behaviors of giant clusters for site percolation 
under random and product rules on two-dimensional lat- 
tice are studied in Sections 3 and 4 respectively. The 
universality classes of bond and site percolation under 
random and product rules are discussed in Section 5. 


II. CRITICAL BEHAVIORS OF PERCOLATION 
PHASE TRANSITION 


The percolation phase transition in a network with 
N = L’ nodes is indicated by the appearance of a giant 
cluster whose size becomes comparable with N. The size 
Sı of the giant cluster is taken as the order-parameter of 
percolation phase transition. Near a critical point, the 
size S2 of the second largest cluster in the network be- 
comes comparable with N also. Both Sı and S2 demon- 
strate critical behavior near the critical point [10, i 

Following the finite-size scaling of order-parameter |23, 
BA, we anticipate that the reduced sizes of S; and Sz ina 
network with N nodes and N, edges follow the finite-size 


scaling form as [10, [1] 


sı (r, L) = S1/L4 = pee 5 (tL/”), (1) 
s2(r, L) = S2/L4 = LP Sa(tL”), (2) 


where r = N,/N is the reduced number of edges and 
t = (r—re)/re is the deviation from critical point re. The 
scaling variable tL*/” is related to the size ratio L/€ of L 
to the correlation length in bulk € = &o|t|~”. The finite- 
size scaling form is supposed to be valid in the asymptotic 
critical region with L > 1 and |t| « 1. 

From Eqs.) and (2), the finite-size scaling form of 
size ratio s2/sı is obtained as 


So/81 = So(tL!”)/51 tL”) =UCLY”). (3) 


At the critical point re, the size ratio s2/sı = U(0) 
and becomes independent of system size L. The curves 
82/8, of different L against r have a cross-point, which 
can be used to determine the critical point re. 

The logarithm of Eq. (I) is 


In sı (r, L) = —(8/v) nL + ng, (tL). (4) 
At the critical point, we have 


In sı(re, L) = —(8/v) nL + ln 51 (0), (5) 


which is a straight line against ln L. We can fix the criti- 
cal point re by the straight line of ln sı (r, L) with respect 
to ln L. From the slope of the straight line, we can de- 
termine the critical exponent ratio 3/v. By introducing 
scaling variable tL!/” to make curves sz /s1 of different 
L collapse into a scaling function U(tL'/”), the critical 
exponent v can be estimated. 

There are critical behaviors also during the evolution 
of a network. The percolation phase transition is accom- 
panied not only by the appearance of a giant cluster but 
also by the giant size jump of the largest cluster. In the 
Refs [17 and Es), the percolation phase transition was 
investigated by the size jumps of the largest cluster. 

In a simulation of network evolution, edges are added 
under some rule one by one into network. At an evo- 
lution step T, the number of edges in the network in- 
creases from T — 1 to T and the size of the largest clus- 
ter varies correspondingly from S(T — 1) to S,(T). At 
this step, the largest cluster has a reduced size jump 
ôr = [Si(T) — S(T — 1)] /N. The largest reduced size 
jump of the whole evolution process is 


A = maz {61, 69, 63,...}. (6) 


The evolution step, where the largest size jump ap- 
pears, is denoted as Te. In the i-th simulation, we obtain 
the largest reduced size jump AČ and the critical re- 
duced evolution step r® = TS? /N. From the results of 
M simulations, we can calculate following averages 


z Pe 
A) = =>, (7) 
i=l 
1 M ; 
Fe(L) = aor. (8) 
j=l 


In the bulk limit L — oo, we suppose that 7,.(L) ap- 
proaches its bulk limit r.(0o) as 


F.(L) = relo0) + ap, L71". (9) 


The bulk limit r.(co) from the network evolution should 
be equal to the critical point re obtained from the largest 
cluster. 

The character of percolation phase transition can be 
determined from the finite-size effect of A(L). If A(L) 
approaches a non-zero value in the limit L — oo, the per- 
colation is a discontinuous phase transition. For a contin- 
uous percolation phase transition, A(L) has a power-law 
finite-size effect as 


A(L) = aa L^. (10) 


In the i-th simulation of network evolution, there are 
fluctuations ôre = r® — 7.(L) and 6A = A® — A(L). 
Their root mean squares are calculated as 

Xr = y< (re)? >, (11) 
xa = V< (6A)? >. (12) 


Their dependence on network size L is described by 
the exponents v and {2 as 


Xr x L71, (13) 
xa x LP, (14) 


In Ref. Ei, the asymptotic behavior of A(L) was used 
to judge the continuity of a percolation phase transition. 
The so-called upper pseudo transition point in Ref. 
was estimated from 7,.(L). The general clique percolation 
phase transition in random networks are identified by 
studying the averages and fluctuations of A and re is). 

We will investigate the above critical behaviors of a 
percolation phase transition, which are characterized by 
the critical exponents 8, 61,82, vV,vı and v2. The differ- 
ent percolation phase transitions with the same critical 
exponents belong to the same universality class. 


II. RANDOM SITE PERCOLATION ON 
TWO-DIMENSIONAL SQUARE LATTICE 


On a two-dimensional square lattice, the percolation 
transition can be investigated both for site and bond. For 
the site percolation, we consider a L x L square lattices 
with periodic boundary conditions in both directions. At 
first, there are N = L x L unoccupied sites in a lattice. 
Then the sites in the lattice are occupied one by one 
under some rule. Clusters consisting of adjacent sites 
appear during the evolution. When N, sites are occupied 
in the lattice, we define a reduced number of occupied 
sites r = N,./N. If the size of the largest cluster becomes 
comparable with N, there is a site percolation in the 
lattice. 

For the random site percolation (SP), sites are added 
randomly. It is known that SP on two-dimensional 
square lattice occurs at re = 0.59276421(13) [19]. From 
Coulomb gas arguments and conformal field the- 
ory [21], the critical exponents 8 = 5/36 and v = 4/3 
were predicted. 

In our Monte Carlo simulation of site percola- 
tion, the algorithm of Newmann and Ziff , 
has been used. We have taken lattice sizes L = 
32, 64, 128, 256, 512,1024,2048,4096 and 8192 in our sim- 
ulations. 12,800,000 independent simulations have been 
run for each lattice size. 

In Fig[I] the size ratio s2/s1 is plotted with respect to 
r for three different lattice sizes. From its cross-point, we 
get the critical point re = 0.59276(3) which agrees with 
the results of Ref. [19}. 

In Fig] Ins; is shown with respect to In L for three 
different r. At r = 0.5927, the curve is nearly a straight 
line with slope equal to —0.108. The curvature becomes 
negative for r = 0.5924 and positive for r = 0.5930. 
We can conclude that the critical point re = 0.5927(3) 
and the critical exponent ratio 6/v = 0.108(4), which 
agrees well with the exact value 8/v = 5/48. Within er- 
ror bounds, the critical point from In sı agrees with that 
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Figure 1: Size ratio of the second largest cluster to the largest 
cluster U = s2/s1. The critical point of site percolation re = 
0.59276(3) is obtained from its cross-point. 
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Figure 2: Log-log plot of sı with respect to L for three differ- 
ent reduced number of occupied sites r. At re = 0.5927(3), 
the curve becomes a straight line and its slope gives the crit- 
ical exponent ratio 8/v = 0.108. 


from s2/s1. Using 1/v = 0.75 in Fig[3] three curves of 
82/8, at L = 256,512, 1024 collapse into a scaling func- 
tion function U(tL/”). 

In Fig] we show critical behaviors of the largest clus- 
ter during the evolution of SP. The Monte Carlo sim- 
ulation results of 7-(L) are demonstrated in Fig/4{A). 
Using re(00) = 0.5927(3) from the above investigation, a 
log-log plot of r-(co) — -(L) with respect to L is made. 
Our Monte Carlo data demonstrate a linear dependence 
of In[r.(co) — re(L)] on nL. From its slope, we obtain 
the inverse of exponent 1/7; = 0.75(1). In Fig/4{B), 
the power-law finite-size behavior of A(L) in Eq.(10) is 
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Figure 3: Finite-size scaling function U(tL'/”) of the size ratio 
82/81 for random site percolation with 1/v = 0.75. 
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Figure 4: Critical behaviors of random site percolation dur- 
ing network evolution: (A)log-log plot of re(o0) — Fe(L) with 
r-(oo) = 0.5927 and the exponent 1/vı = 0.75(1); (B)log-log 
plot of A(L) with the exponent 82 = 0.104(1); (C)log-log plot 
of Xr with the exponent 1/v2 = 0.74(1); (D)log-log plot of ya 
with the exponent $2 = 0.104(1). 


confirmed by the Monte Carlo simulation data. The ex- 
ponent 3; = 0.104(1) is calculated from the slope of the 
straight line. 

The root mean squares of fluctuations y, and ya are 
shown in Figs[4{C) and (D). Our Monte Carlo simulation 
data of Xr follow a power-law dependence on L with the 
exponent 1/v2 = 0.74(1). The power-law behavior of ya 
in Eqff4]is confirm by our simulation data and we obtain 
Bq = 0.104(1). 

We have investigated also the evolution critical be- 
haviors of random bond percolation (BP) on two- 
dimensional square lattice. These results of BP are given 
in Tablefl] We can see that the critical exponents of SP 
are equal to that of BP. Further, the critical exponent 
ratio 8/v is equal to 6; and 82 and v is equal to vı and 
və within error bounds. 


IV. SITE PERCOLATION UNDER PRODUCT 
RULE ON TWO-DIMENSIONAL SQUARE 
LATTICE 


In this section, we consider the site percolation under 
product rule (SPPR). We take the product rule [14, 
for the site evolution as follows: (1)selecting two unoccu- 
pied sites randomly; (2)calculating the product of sizes 
of the clusters which are connected by two chosen sites 
respectively; (3)the site with smaller product is occupied. 

The reduced sizes sı of the largest cluster for lattice 
size L from 512 to 8192 are shown in Fig[5} sı is nearly 
zero at small r and becomes finite when r is larger than 
a critical value re. Finite sı indicates the emergence of a 
giant cluster and a site percolation transition. 

The results obtained from the size of the largest clus- 
ter are shown in Figs. [6] [7] and B| From the fixed 
point of s/s, in Fig] we can get the critical point 
re = 0.7634(2). The straight line in Fig[7] gives the crit- 
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Figure 5: Reduced size of the largest cluster as a function of 
r at lattice size L = 512, 1024, 2048, 4096, 8192. 
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Figure 6: Size ratio of the second largest cluster to the largest 
cluster U = s2/si at different lattice sizes. Its fixed point 
gives the critical point re = 0.7634(2). 


ical point re = 0.7634(5), which is in agreement with 
the value of the fixed point. The accuracy of our critical 
point is higher than that of re = 0.768(3) in Ref.[14] and 
re = 0.756(6) in Ref. [15]. 

The slope of the straight line in Fig[/]gives the critical 
exponent ratio 6/v = 0.010(1), which is different from 
B/v = 0.04 + 0.02 of Ref. [15]. 

We can see in Fig[8] that the size ratios s2/sı of dif- 
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Figure 7: Ins; with respect to In L around the critical point 
re = 0.7634(5), where the curve becomes a straight line with 
slope equal to —0.010. 
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Figure 8: Finite-size scaling function U(tL'/”) of size ratio 
82/81 with 1/v = 0.42. 


ferent lattice sizes collapse well with the scaling variable 
tL'/” for data of L > 512. This means that the asymp- 
totic region of SPPR is more restricted than that of SP. 
In the finite-size scaling plot of s2/s1, we use the inverse 
of the critical exponent 1/v = 0.42(1). 


The size ratio at critical point U(0) is supposed to be 
universal. We obtain U(0) = 0.03 for SPPR and U(0) = 
0.1 for SP. This difference of U(0) indicates also that 
SPPR and SP belong to different universality class. 


We have studied also the size jumps of the largest 
cluster during network evolution for SPPR. The differ- 
ent quantities related to the size jump are shown in 
Fig[9] with respect to the lattice size L. From Fig[9{A), 
we can conclude that the average transition point ap- 
proaches to its bulk limit according the power-law in 
Eqfii] Using the Monte Carlo data of the lattice size 
l > 512, we get the inverse of exponent 1/1; = 0.42(2). 
The root mean square of the fluctuations of transition 
point is plotted in Fig[{C) and a power-law in Eq[13]is 
found 1/12 = 0.42(3). The simulation results of the size 
jumps of the largest cluster are presented in Figs[9{B) 
and (D). The average size jump A shows a quite differ- 
ent L-dependence from that of the random site percola- 
tion. No power-law behavior is found for A(L), which 
increases with L at first and then decreases. The root 
mean square of fluctuations of A follows a power-law of 
the lattice size also. From Monte Carlo simulation data, 
we get 32 = 0.008(4). 


For comparison, we have studied also the evolution 
critical behaviors of bond percolation under product rule 
(BPPR). The critical exponents related to the cluster 
size jump 6, = 0.044(5) and 2 = 0.041(4) are ob- 
tained. For the exponents related to critical point, we 
get 1/vı = 0.99(5) and 1/v2 = 0.94(3). The critical expo- 
nent ratio 8/v and the inverse of critical exponent 1/v of 
BPPR have been calculated in Ref.[10). 8/v = 0.064(3) 
is different from the exponents 8; and 62. Within error 
bounds, 1/v = 0.93(1) is equal to 1/vı and 1/12. All 
results of BP, SP, BPPR and SPPR are summarized in 
Table[I] We can see that the critical exponents of BPPR 
are different from that of SPPR. The universality classes 
of BPPR and SPPR are different. 
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Figure 9: Results of site percolation under the product rule: 
(A)log-log plot of r-(co) — 7-(L) with re(co) = 0.7634 and 
the exponent 1/1, = 0.42(2); (B)log-log plot of A(L), which 
increases at first and then decreases slowly; (C)log-log plot of 
Xr with the exponent 1/v2 = 0.42(3); (D)log-log plot of ya 
with the exponent 62 = 0.008(4). 


Table I: Summary of critical points and critical exponents of 
BP, SP, BPPR and SPPR. re, 8/v and 1/v of BPPR are taken 
from Ref. [10]. 


BP SP BPPR SPPR 
re 0.5000(4) 0.5927(3) 0.5266(1) 0.7634(2) 
1/v 0.75 0.75(1) 0.93(1) 0.42(1) 
1/1 0.75(1) 0.75(1) 0.99(5) 0.42(2) 
1/v2 0.74(1) 0.74(1) 0.94(3) 0.42(3) 
B/v 0.108(4) 0.108(4) 0.064(3) 0.010(1) 
Bı 0.105(2) 0.104(1) 0.044(5) 
B2 0.104(1) 0.104(1) 0.041(4) 0.008(4) 


V. CONCLUSION 


We have studied the critical behaviors of the sizes and 
the size jumps of giant clusters on two-dimensional square 
lattice for both site and bond percolation under ran- 
dom and product rules respectively. Under the random 
rule, sites or bonds are added randomly into the lattice. 
We add the site or bond with smaller size product of 
two connected clusters into the lattice from two candi- 
dates picked randomly when the product rule is taken. 
The finite-size scaling forms of giant cluster size and size 
jumps are proposed and used to determine the transition 
points and critical exponents of the percolation transi- 
tions from Monte Carlo data. 

Our results show that the critical exponents of the size 
and size jump of random site percolation are equal to 
that of random bond percolation within error bounds of 
Monte Carlo data. As expected, the universality class 
of percolation is independent of site and bond under the 
random rule. The critical exponents of the size and size 
jump in the site percolation under the product rule are 
different from that in the random percolation and also 
the bond percolation under the product rule. The site 


percolation and bond percolation under the product rule 
do not belong to the same universality class. 
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